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ABSTRACT 

Let K be a field and (/i, . . . , f n ) C K[Ai, . . . , X n ] be a 
sequence of quasi-homogeneous polynomials of respective 
weighted degrees (di, . . . ,d n ). By this, we mean that there 
exists (wi, . . . , w n ) £ Z™o s.t. for any 1 < j < n, the poly- 
nomial /^(XJ" 1 , . . . , X™ n ) is homogeneous and has degree 
di. This paper focuses on designing strategies for comput- 
ing Grobner bases of such systems. 

Under some genericity assumptions on the coefficients of 
the fi's, we prove the bound (X^ILiC^* ~ w i)) +max{u>i} on 
the degree of regularity of such systems. This is done by us- 
ing well-known results on Hilbert series and designing a ded- 
icated quasi-homogeneous matrix-Fs algorithm taking care 
of this structure. We provide a detailed complexity analysis 
of this algorithm, and we prove that we can divide the com- 
plexity bounds by P u , where P — Y\2=i Wi - We also show 
how to use the FGLM (change of ordering) step in this spe- 
cial setting to obtain a Grobner basis for the lexicographical 
ordering in time polynomial in the weighted Bezout bound 

nn 
1=1 ~wl' 

All in all, we provide a complete computational strategy 
for solving quasi-homogeneous systems in time polynomial 
in this weighted Bezout bound. We finally prove that these 
complexity results apply for O-dimensional affine systems 
which are obtained by specializing to 1 a variable in a quasi- 
homogeneous system. We provide some experiments show- 
ing that this computational strategy sometimes offer gains 
of up to 10, and allows to handle systems which cannot be 
tackled with previous approaches. 

1. INTRODUCTION 

Motivation and problem statement. Polynomial systems 
solving is a very important problem in computer algebra, 
with a wide range of applications in theory (algorithmic ge- 
ometry) or in real life (cryptography). For this purpose, 
Grobner bases of polynomial ideals are a valuable tool, and 
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practicable computation of the Grobner bases of any given 
ideal is a major challenge of modern computer algebra. Since 
their introduction in 1965, many algorithms have been de- 
signed to compute Grobner bases ([6, 9, 10, 13]), allowing 
for more and more systems to be solved. 

Systems arising from "real life" problems often show some 
structure easing the Grobner basis computation. For ex- 
ample, it is known that homogeneous systems, or system 
with an important maximal homogeneous component, are 
better solved by using a degree-compatible order, and then 
applying a change of ordering. In this paper, we study 
a structure slightly more general than homogeneity, called 
quasi-homogeneity. More precisely, we will say that a poly- 
nomial P(Xi, . . . ,X n ) is quasi-homogeneous for the system 
of weights W = [wi, . . . , w n ), if the polynomial 

Q(Yi,...,Y n ) := P(Y?\...,Y?~) 
is homogeneous. Systems with such a structure are likely 
to arise for example from physics, where all measures are 
associated with a dimension which, to some extent, we can 
see as a weight. 

Let F = (fx, ...,/m) be a zero-dimensional system of 
quasi-homogeneous polynomials (or of polynomials with an 
important maximal quasi-homogeneous component), it is 
possible to try to compute directly a Grobner basis of the 
ideal generated by F, ignoring the quasi-homogeneous struc- 
ture. However, there is no general way of evaluating the 
complexity of this strategy. 

Another approach is to compute the homogenized system 
F :— {fiiX™ 1 , . . . , X™"))and then compute a Grobner basis 
of that system, using the usual strategies for the homoge- 
neous structure. The first step of the computation proves 
much faster than with the naive strategy, but the change of 
ordering is slowed down, and thus becomes the main bottle- 
neck of the computation. This leads to the following three 
questions : 

1. Is this strategy valid: does a Grobner basis of the ideal 
generated by F give any information regarding a Grob- 
ner basis of {F}7 

2. Can we obtain a similar speedup for the change of 
ordering? 

3. Can we evaluate the complexity of this strategy, and 
thus explain these speedups? 

Main results. For the strategy consisting in ignoring the 
quasi- homogeneous structure of the system, there is no pre- 
cise complexity estimate. On the other hand, in this paper, 



we offer a full strategy to compute a lexicographical Grobner 
basis of a quasi-homogeneous zero-dimensional system, with 
complexity results under genericity assumptions, so that the 
answer to the three questions above is positive. 

More precisely, we define in a natural way a weighted 
variant of the GRevLex ordering, which we call the W- 
GRevLex ordering. We then prove that a GRevLex Grob- 
ner basis of (F) provides a VF-GRevLex Grobner basis of 
(F), and we can use the FGLM algorithm to perform the 
change of ordering to the Lex ordering, without any slow- 
down. For complexity results, we will assume the system to 
satisfy the two following generic conditions: 
HI. The sequence fi, ■ ■ ■ , fm is regular; 
H2. The sequence is in Noether position w.r.t. 

Xi, . . . , Xi, for any 1 < i < m. 

Under the hypothesis HI, we adapt the known results 
of the homogeneous case, using similar arguments based on 
Hilbert series, to give estimates on the degree of the ideal 
and the degree of regularity of the system : 

deg(J) = Y\ — ; d Icg (F) < ^2 (di ~ w ») + max{w;}. 

i=l W ' i=l 

We also present a variant of the matrix-Fs algorithm which 
follows more closely the behavior of F5 ran on F. A combina- 
torial result found in [I] shows that the number of columns 
of the matrices appearing in that variant of matrix-Fs is 
approximately divided by Yl7=i Wi > wnen compared to the 
regular matrix-Fs algorithm. In the end, we show that 
for systems satisfying HI, our strategy, using our variant 
of matrix-Fs and FGLM, performs in time polynomial in 
n™=i di/ nr=i Wi > that is polynomial in the number of solu- 
tions. 

Further assuming the hypothesis H2, we also carry on the 
precise complexity analyses done in [2] for homogeneous sys- 
tems, and adapt them to the quasi-homogeneous case to de- 
duce a precise complexity bound for our quasi-homogeneous 
variant of Matrix-Fs ■ It confirms that overall, we are able to 
compute a Lex Grobner basis for a generic quasi-homogeneous( 
system P w times faster than for a generic homogeneous sys- 
tem, where P = Yl"=i w i- 

We have been able to run benchmarks on both generic 
systems and real-life systems arising in cryptography. Ex- 
perimentally, in both cases, our strategy seems always faster 
than ignoring the quasi-homogeneous structure, the gain in- 
creasing with the considered weights. 

Experiments have also shown that the order of the vari- 
ables can have an impact on the performances of both strate- 
gies, predicting it seems to require more sophisticated tools 
and is left to further work. 

Prior works. Making use of the structure of polynomial 
systems to develop faster algorithms has been a general 
trend over the past few years: see for example [14], [7] or 
[15]. Polynomial algebras graded with respect to a system 
of weights have been studied from the sight of commutative 
algebra. Most notably, the Hilbert series of ideals defined 
by regular sequences, which we use several times in this pa- 
per, is well known, and could be found for example in [19]. 
The paper [ 1 8] defines many structures of polynomial alge- 
bras, including weighted gradings, in preparation for future 
algorithmic developements. Combinatorial objects arising 
when we try to estimate the number of monomials of a given 
W-degree are called Sylvester denumerants, and studied for 
example in [1]. 



When it comes to Grobner bases, weighted gradings and 
related orderings have been described in early works such 
as [4]. However, as far as we know, the impact of a quasi- 
homogeneous structure on the complexity of Grobner bases 
computations had never been studied. 

Amongst the various computer algebra software able to 
compute Grobner bases, it seems only Magma [5] is taking 
advantage of a system of weights, and if available, uses the 
V^-GRevLex ordering to compute an intermediate basis be- 
fore the change of ordering. Moreover, it is able to retrieve 
the appropriate system of weights. However, this strategy is 
only available for quasi-homogeneous systems. 
Structure of the paper. In section 2 we define more pre- 
cisely quasi-homogeneous systems, and we compute their de- 
gree and degree of regularity assuming the above hypotheses. 
We also take this occasion to quickly show that these hy- 
potheses are generic. In section 3 we prove that the strategy 
consisting of modifying the system is correct, we explain how 
we can adapt matrix-Fs and FGLM to quasi-homogeneous 
systems, and then we evaluate the complexity of these algo- 
rithms. In section 4, we briefly explain how these results for 
quasi-homogeneous systems still help in case the system is 
obtained from a quasi-homogeneous system, by specializing 
one of the variables to 1. We also give an example of such 
a structure, as well as the associated algorithm. Finally, in 
section 5, we give some experimental results. 

2. QUASI-HOMOGENEOUS SYSTEMS 
2.1 Weighted degrees and polynomials 

Let K be a field, we consider the algebra A := KLYi , . . . , X n ] 
= K[X]. Even though usually, one uses the total degree to 
grade the algebra A, there are other ways to define such a 
grading, as seen in [1], for example. 

Definition 1. Let W = (wi, . . . ,w n ) be a vector of posi- 
tive integers. Let a = (ai, . . . , a n ) be a n-uple of nonnega- 
tive integers, we call W -degree, or weighted degree of a mono- 
mial X Q = A" 1 • • • X% n the integer deg w (X Q ) = £™ =1 Wi m!g 
The vector W is called system of weights. We denote by 1 
the weight system defined by (1, ... , 1), associated with the 
usual grading on A. 

One can prove that any grading on K[X] comes from such 
a weight system ([4, sec. 10.2]). We denote by (K[X],W) 
the VF-graded structure on A, and in that case, to avoid am- 
biguities, we use the adjective W -homogeneous for elements 
or ideals, or quasi-homogeneous if W is clear in the context. 
The word homogeneous will be reserved to 1-homogeneous 
items. 

Proposition 1. Let (K[Xi, . . . ,X„], W) be a graded poly- 
nomial algebra. Then the application 

hom w : (K[Xi,...,X n ],W) -)■ (K[ti,. ..,«„], 1) 
/ ^ /(C, ■■■>*»") 

is an injective graded morphism, and in particular the image 
of a quasi- homogeneous polynomial is a homogeneous poly- 
nomial. 

Proof. It is an easy consequence of the definition of the 
weighted grading. □ 

The above morphism also provides a quasi-homogeneous 
variant of the GRevLex ordering, which we call the W- 
GRevLex ordering: 

u <iv- g rcvi cx v hom w (v) < grov icx hom w (v) 



Given a W-homogeneous system F, one can build the ho- 
mogeneous system hom w (F), and then apply classical al- 
gorithms ([10, 13]) to that system to compute a GRevLex 
(resp. Lex) Grobner basis of the ideal generated by hom lv (i ? )| 
We will prove in section 3 (prop. 7) that this basis is con- 
tained in the image of hom w , and that its pullback is a 
VF-GRevLex (resp. Lex) Grobner basis of the ideal gener- 
ated by F. 

Let us end this paragraph with some notations and def- 
initions. We call degree of regularity of the system F the 
highest degree d Icg (F) reached in a run of F5 to compute 
a GRevLex Grobner basis of hom w (F). We call index of 
regularity of an ideal I the degree i rog of the Hilbert series 
HSa//, defined as the difference of the degree of its numer- 
ator and the degree of its denominator. 

Recall that given a homogeneous ideal I, we define its 
degree D as the degree of the projective variety V(I), as 
introduced for example in [16]. This definition still holds for 
the quasi- homogeneous case. In case the projective variety 
is empty, that is if the affine variety is equal to {0}, we 
extend that definition by letting D be the multiplicity of 
the point, that is the dimension of the K-vector space A/F 
Finally, from now on we will only consider affine varieties, 
even when the ideal will be quasi-homogeneous. It means in 
particular that the dimension of V(0) is n, and that a zero- 
dimensional variety will be defined by at least n polynomials. 

2.2 Degree and degree of regularity 

Zero-dimensional regular sequences. As in the homoge- 
neous case, regular sequences are an important case to study, 
because it is a generic property which allows us to compute 
several key parameters and good complexity bounds. We 
first show how we can compute the degree, and a bound on 
the degree of regularity of a zero-dimensional ideal defined 
by a regular sequence. 

Theorem 2. Let W — (wi, . . . , w n ) be a weight system, 
and F = (/1, . . . , f m ) a regular sequence of W -homogeneous 
polynomials, of respective W -degrees di,... ,d m . Further as- 
sume that the set of solutions is zero-dimensional, that is 
m = n. We denote by I the quasi-homogeneous ideal gener- 



Furthermore, the degree of the ideal / is equal to the di- 
mension of the vector space A/ 1, that is the value of the 
Hilbert series at t = 1, that is FT" , — ■ □ 

Note that except for this inequality, not much is known 
about the degree of regularity of a quasi-homogeneous sys- 
tem. In particular, the above bound is nothing more than 
a bound, even in the generic case. Let me introduce some 
examples of the three cases one can observe with a quasi- 
homogeneous generic system: 

1. W = (3, 2, 1), generic system of VF-degree D = (6, 6, 6): 



then d T , 



+ 1 = 13; 



2. W = (1, 2, 3), generic system of VF-degree D = (6, 6, 6): 
then dreg = 15 > i rcg + 1 = 13; 

3. W = (2,3), generic system of W-degree D = (6,6): 



then gL 



6 < i r , 



7. 



ated by F. Then we have deg(J) = Yll 



andd Icg (F) <£™=i 



Proof. We will read the degree and degree of regularity 
of the system on the Hilbert series (or Poincare series) of 
the algebra A/F A classical result which can be found for 
example in [19, cor. 3.3] states that, for regular sequences, 
this series is 

/n ; (l-t<"i)---(l-i™«) v ' 

We assumed n — m, so the Hilbert series can be rewritten 



ns A/1 (t) 



(1 + 



dl-lN 



(! + • 



+ t d - 1 ) 



(H h^i- 1 ) . . . (H h t 1 ""- 1 ) ' 

In the 0-dimensional case, recall that the Hilbert series is 
actually a polynomial, and has degree i rog = $2™ =1 (di — Wi). 
This means that all monomials of VF-degree greater than 
i reg are in the ideal, and as such, that the leading terms of 
the VF-GRevLex Grobner basis of F need to divide all the 
monomials of W-degree greater than i rog .This proves that all 
the polynomials in the Grobner basis computed by F5 have 
W-degree at most i reg +max{wi}. And since the F5 criterion 
([10]) ensures that there is no degree-fall in a run of F5 on 
a regular sequence, the algorithm indeed stops in degree at 
most i reg + m&x{wi}. 



Only the case 1 is observed with generic homogeneous sys- 
tems. Furthermore, examples 1 and 2 show that the degree 
of regularity depends on the order of the variables (chosen 
in the description of the weight system). As the Hilbert se- 
ries of a generic sequence doesn't depend on that order, it 
shows that we probably need to find a better tool in order to 
evaluate more precisely the degree of regularity in the quasi- 
homogeneous case. However, the above bound already leads 
to good improvements on the complexity bounds, as we will 
see in the following sections. Also note that these compu- 
tations only hold when the system is 0-dimensional, we will 
discuss that restriction in section 2.3. 

Genericity. We now prove that zero-dimensional W-homo- 
geneous sequences of given W-degree are generically regular, 
under some assumptions on the W-degree. Let's start with 
the first part of this statement: 

Lemma 3. Let n be a positive integer, and consider the 
algebra A := K.[Xi, . . . ,X n ], graded with respect to the sys- 
tem of weights W = (wi, . . . ,w„). Regular sequences of 
length n form a Zariski-open subset of all sequences of quasi- 
homogeneous polynomials of given W-degree in A. 

Proof. Let (di,...,d m ) be a family of W-degrees, we 
consider the set V (K[a] [X]) of all systems of quasi-homogeneous] 
((feeljo^^naisc^fJ^ degree di , . . . , d m . We denote by /1 , . . . , / m | 

the polynomials 01 the generic system, and by I the ideal 
they generate, in K[a][X]. 

Since the Hilbert series (1) characterizes regular sequences 
([19, cor. 3.2]), the sequence (/;) is regular if and only if 
the ideal I contains all monomials of W-degree between 
i Yes {I) + 1 and i les (I) + maxjwi}, where i reg (/) is given 
by ~^2(di — Wi). This expresses that a given set of linear 
equations have solutions, and so it can be coded as some 
determinants being non-zero. □ 

There are some systems of W-degree for which there is 
no regular sequence. The reason is that because of the 
weights, for some systems of IF-degrees, there exists no or 
very few monomials. For example, take n = 2, W — (1,2) 
and D = (1,1). All quasi- homogeneous polynomials of W- 
degree 1 are in M.X, so there is no regular sequence of quasi- 
homogeneous polynomials with these W-degrees. 

However, if we only consider "reasonable" systems of W- 
degrees, that is systems of IF-degrees for which there exists 
a regular sequence, regular sequences form a Zariski-dense 
subset from the above. 

Remark 1. A sufficient condition for example is to take W- 
degrees such that Wi is divisible by d\ , . . . , w n is divisible by 



d n . This way we can define the sequence 

which is regular, and so for such systems of weight, the reg- 
ularity condition is generic. 

We only proved the genericity for quasi-homogeneous se- 
quences of length n, the more general case of a sequence of 
length m < n will be proved in section 2.3 (remark 2). 

2.3 Noether position 

To compute the degree and degree of regularity of quasi- 
homogeneous systems of positive dimension, we will assume 
that the system F = . . . , f m ) we consider is in Noether 
position (as seen in [8, ch. 13, sec. 1] or [3, def. 2]), i.e. 

• for i < m, the canonical image of Xi in K[X]/J is an 
algebraic integer over K[X m +i , . . . , X„]; 

• K[x m +i,...,x„]ni = 0. 

Lemma 4. Let F = /i, . . . , f m be a regular quasi-homo- 
geneous sequence of polynomials in K.[X%, . . . ,X n ]. The se- 
quence F is in Noether position if and only if the sequence 
-Foxt := fi, ■ ■ - , fm, X m+1 ,...,X n is regular. 

Proof. Let I be the ideal generated by the /;'s. The ge- 
ometric characterization of Noether position (see e.g. [17]) 
shows that the canonical projection onto the m first coordi- 
nates 7T : V(I) — > V({Xi, . . . , Xm)) is a surjective morphism 
with finite fibers. This implies that the variety V({F ext )), 
that is 7r _1 (0), is zero-dimensional, and so the sequence is 
regular. 

Conversely, assume _F cx t is a regular sequence. We want 
to show that every Xi for i < m is integral over the ring 
KLY m +i, . . . , X n ]. Since it defines a zero-dimensional ideal, 
for any i, there exists ni £ N such that X™* = LT(/) with 
/ £ (fcxt) for the GRevLex ordering with X\ > ■ ■ • > X„. 
By definition of the GRevLex ordering, we can assume that 
/ simply belongs to /. This shows that every Xi is integral 
over K[Xi+i, . . . ,X n ]/L We get the requested result by in- 
duction on i : first, this is clear if i = m. Now assume that 
we know that K[Xi, . . . , X n ]/I is an integral extension of 
KLY m +i, . . . , X n ]- From the above, we also know that 
is integral over K[Xi, . . . , X n ], and so, since the composition 
of integral homomorphisms is integral, we get the requested 
result. 

Finally, we want to check the second part of the defini- 
tion of Noether position. Assume that there is a non-zero 
polynomial in K[X m+ i , . . . , X n ] H /, since the ideal is quasi- 
homogeneous, we can assume this polynomial to be quasi- 
homogeneous. Either this polynomial is of degree 0, either 
it is a non-trivial syzygy between X m +i, . . . ,X n , so either 
way, it contradicts the regularity hypothesis. □ 

As we did for regular sequences, we first show how we can 
evaluate the degree and degree of regularity of a sequence 
in Noether position, and then we show that the Noether 
position property is generic under some assumptions on the 
W-degree of the polynomials. 

Theorem 5. Let W = (wi, . . - , w n ) be a weight system, 
and fi , . . . , f m a regular sequence in Noether position, of 
quasi-homogeneous polynomials of W -degrees (di, . . . , dm)- 
The same way we did above, we denote by I the ideal gen- 
erated by the fi 's. Then we have deg(J) = YliLi ■ an< ^ 
d reg (7) < YXLi i d i ~ w i) + max{u>i}. 

Proof. Let us denote by I' the ideal generated by -F ex t- 
The degree of the ideal /' is the same as that of I, because 
the variety it defines is the intersection of V(I) with some 



non-zero-divisor hyperplanes. Furthermore, all critical pairs 
appearing in a run of F5 on F will also appear in a run of 
F 5 on F cxt , ensuring that d rcg (F) < d rcg (F cxt ). 

But since by Noether position, the family F c ^t defines a 
zero-dimensional variety, we can use the previous computa- 
tions to deduce its degree of regularity and the degree of 
/'. □ 

Lemma 6. Let n be a positive integer, and consider the 
algebra A := K[Xi, . . . ,X n ], graded with respect to the sys- 
tem of weights W = (wi, . . . , w„). Systems in Noether po- 
sition form a Zariski-open subset of all systems of quasi- 
homogeneous polynomials of given W -degrees in A. 

Proof. Let F — . . . , fm) be m generic quasi-homoge- 
neous polynomials, with coefficients in K[a]. We use the 
same characterization of a zero-dimensional regular sequence 
as we did in the proof of Lemma 3. It allows us to express the 
regularity condition for the sequence (/1, . . . , f m , X m +i, ■ ■ ■ , 
as some determinants being non-zero, which by definition, 
shows that the condition of being in Noether position is an 
open condition. □ 

Since a sequence in Noether position is in particular a reg- 
ular sequence, we are confronted with the same problem as 
for the genericity of regular sequences, that is the possible 
emptiness of the condition. However, it is still true that for 
"reasonable" systems of W-degrees, that is systems of W- 
degrees for which there exists enough monomials, sequences 
in Noether position do exist, and thus form a Zariski-dense 
subset of all sequences. For example, since the sequence 
Ij 1 '™ 1 ,.--^™"'"'" is m Noether position, the sufficient 
condition given in Remark 1 is also sufficient to ensure that 
sequences in Noether position are Zariski-dense. 

Remark 2. Any sequence in Noether position is in partic- 
ular a regular sequence, so Lemma 6 proves that, under the 
same assumption on the degree, regular sequences of length 
m < n are generic amongst quasi-homogeneous sequences of 
given V^-degree. 

3. COMPUTING GROBNER BASES 
3.1 Using the standard algorithms on the ho- 
mogenized system 

As we said before, in order to apply the F5 algorithm to 
a quasi- homogeneous system, we may run it through hom w . 
This is shown by the following proposition. 

Proposition 7. Let F = (/1, . . . , f m ) be a family of poly- 
nomials in K[Xi, . . . , X n ], assumed to be quasi-homogeneous 
for a weight system W = (wi, . . . , w n ). Let <i be a mono- 
mial order, G the reduced Grobner basis of hom w (F) for 
this order, and <2 the pullback of <i through hom^,. Then 

1. all elements of G are in the image of hom w ; 

2. the family G' := homjT / 1 (G) is a reduced Grobner basis 
of the system F for the order <2. 

Proof. For the first statement, it is enough to notice that 
we can compute such a reduced Grobner basis by applying 
the Buchberger algorithm and reducing that base. Both 
steps only involve computations of S-polynomials, and one 
checks that S-Pol(hom w (/), hom w (g)) — hom M/ (S-Pol(/, 
So at the end of the computations, all polynomials in the ba- 
sis are in the image of hom w . 

For the second assertion, let us first notice that since <2 is 
the reverse image of <i by hom w , one has, for any polyno- 
mial /, homr^(LT(/)) = LT(hom7 v 1 (/)) and so G' is indeed 



a Grobner basis of (G 1 ). What remains to check now is that 
this ideal is equal to (F). This makes use of the following 
lemma: 

Lemma 8. Let P — (pi, . . . ,p^) be a family of polynomi- 
als, then hom^ 1 ((hom w (P))) = (P). 

Proof. The reverse inclusion is easy: if/ = ^2liPi £ (P)( 
then hom w (f) = ]T hom w (qi) hom^p,) € (hom w (P)). 
Let's prove the direct inclusion. Let g be a polynomial in 
hom^, 1 ((hom w (P))) , there exist polynomials gt such that 
hom ff ( 9 ) = E9.hom w .(p,). 

Separating the monomials between those which are in 
hom w/ (K[X]) and those which aren't, we have a decompo- 
sition of all gi's as gi — hom w (qi) + r\, where no mono- 
mial of n is in hom M/ (K[X]). So we have the decomposition 
hom w (p) = hom w (Y.qiPt) + J2 nh.om. w (j>i). 

So we want to show that the second term in that sum is 
zero. This is true, because if two polynomials p and q are 
made only of monomials out of hom lv (K[X]), their sum is 
either also made only of such monomials, either zero. So here 
we have rihom w (pi) = 0, and so hom w (g) = hom w (£] q%Pi, 
Since the morphism hom w is injective, this shows that g — Y <li 

Applying that lemma twice, one checks that 
<G') = hom^ 1 ({hom w (G'))) = hom w 1 ({G)) 

= hom w . 1 ((hom M ,(F))) = (F), 
which is the wanted result. □ 

In practice, if we want to compute a Lex Grobner basis 
of F, we generate the system F — hom w (F), we compute a 
GRevLex basis Gi of F with F5 or matrix-Fs, and then we 
compute a Lex Grobner basis G2 of F with FGLM. In the 
end, we get a Lex Grobner basis of F, which we turn into a 
Lex Grobner basis of F via hom^ . 

We immediately notice that this method has several draw- 
backs for both matrix-Fs and FGLM: we have seen, at least 
for regular sequences, that going through hom^ increases 
the degree and the known bound for the degree of regularity 
of the system. Furthermore, for a given integer d, there are 
way more monomials of 1-degree d than of W-degree d. This 
results in matrix-Fs working with matrices that are larger 
than necessary. 

3.2 Direct algorithms 

Let us first consider the problem of the change of ordering. 
In the above process, we can apply hom^, 1 to the basis Gi 
and this way obtain a W-GRevLex basis Gi of F. We can 
then run FGLM on that basis to obtain a Lex basis of F. 
This way, we can avoid the problem of a greater degree of 
the ideal on the complexity of the FGLM step. 

About matrix-Fs, we can also avoid going through hom w , 
all we have to do is to change the algorithm a little, in 
order to directly consider the variables with their weight. 
The modified algorithm is algorithm 1 opposite. The func- 
tion F5Criterion(^i, i, M) implements the F 5 -criterion de- 
scribed in [10]: it evaluates to false if and only if fi is the 
leading term of a line of the matrix Md-d it i-i- The function 
EchelonForm(M) reduces the matrix M to row-echelon 
form, not allowing any row swap. 

The only part that differs from the classical matrix-Fs is 
the loop, from line 9 to line 15. Recall that the classical al- 
gorithm loops on the labels (e, fi) in degree d— 1, and check 
the F5 criterion on all new monomials of degree (d — deg(/j)) 
one can obtain from e. Instead, this algorithm loops on 



Algorithm 1: Matrix-F 5 (W-homogeneous version) 



Input: 



3 
4 
5 

6 
7 

8 
9 

10 
11 
12 
.13 



I 

Phi 



fx , . . . , f m W-homogeneous polynomials 

^max £ N 

Output: G Grobner basis of . . . , f m ) up to 
W-degree d max 

1 <?«- ; 

2 for d = 1 to d max do 
Md,Q matrix with lines; 
for i = 1 to m do 

if d — deg(/i) then 

I M d , t «- M d ,i-i U ligne fi with label (1, f z ); 
else if d > deg(/i) then 
M d ,t <- Md,i-i ; 

foreach n monomial of VF-degree d — deg(/;) 
do 

if F5Criterion(/^, i, M) then 

Xj 4— biggest variable dividing /1 ; 

e <- fi/xj ; 

if there exists a line of label (e, fi) in 
Md- Wj ,i then 

/ reduced polynomial 
associated to that line; 
Md,i Md,i U Xjf with label 

re Md.m «- EchelonForm (Md,m ) ; 
17 For any line having been reduced to a non-zero 

polynomial, append it to G ; 
is return G 



monomials of W-degree (d — deg(/i)), checks the F5 crite- 
rion on these monomials, and then only checks if there's a 
label (e, fi) (possibly in W-degree lower than d — 1) which 
might lead to that monomial. 

Note that the F5 algorithm with critical pairs selects the 
polynomials it reduces in a finer way than the matrix-Fs al- 
gorithm. In the W-homogeneous case, applied on F, it will 
automatically "jump" from W-degree d to W-degree d + Wi, 
without considering the parasite monomials. This gives an 
explanation to the fact that the F5 algorithm seems to run 
much faster on homogeneous systems obtained through hom^l 
than on generic homogeneous systems. 

3.3 First complexity bounds 

Let F = (/1, . . - , fn) be a system of W-homogeneous poly- 
nomials in K[X\ , . . . , X n ] , and let I be the ideal generated 
by F, D the degree of I, d rog its degree of regularity and i reg 
its index of regularity. The classical complexity bounds of 
matrix-Fs (for a regular system) and FGLM are 

Cf 5 = O {d ICS Md TBIl ,w(nT) ; Cfglm = O (nP> 3 ) , 
where Md,w(n) stands for the number of monomials of W- 
degree d in n variables. 

Assuming the system F is a regular sequence, we have 
already seen the following estimates:d lcg < i rcg + max{ii)i} 
and D = Yl7=i &il XS2=\ Wi - ^ we compare these values with 
their equivalent with the system of weights 1, we notice an 
important gain in complexity in the FGLM algorithm, as 
well as a significant practical gain in matrix-F 5 . 

But the real gain in matrix-Fs is way more important, be- 
cause the matrices of polynomials of given W-degree are 
much smaller than the matrices of polynomials of given 



1-degree. The point of the following lemma is to evaluate 
this difference. 

Lemma 9. Let W = (wi, . . . , w n ) be a weight system, and 
for any i, let Wi = (wi, . . . ,Wi). For any integer d, we de- 
note by Md,w(n) the number of monomials of W -degree d, 
that is the size of the matrix of W -degree d. Let 8 := gcd (W), 
P := n"=i Wi > the integer defined recursively as following: 

a n a a , gcd (Wj-i) 

Si = 0, Si = Si-i + Wi ■ , lxr for i>2 

gcd (Wi) 

and Ti the integer defined recursively as following: 



Ti =0, Ti = Ti^i+Wi- 



;cd (Wi- 



M - 1 ) - 1 for i > 2. 



gpd(Wi) 

Then the number of monomials of W -degree d is bounded 
above and below by: 

s s 

— M d -T n - n +i,x(n) < M d ,w(n) < — M d +s n -n+i,x(n). 

Proof. This is a consequence of theorems 3.3 and 3.4 
in [ ], if we recall that M d ,x(n) = ( d+n d ~ 1 ) = C^ 1 )- □ 

Note that if W = 1, the bounds we get are trivial, which 
means the complexity bounds we will obtain with them will 
specialize without any difficulty to the known bounds for the 
homogeneous case. 

Using the notation S — S n , we get this new complexity 
bound for quasi-homogeneous matrix-Fs: 
CV 5 = O {d ICS M dl ^ w {nY) 



O 



{m}) 



max 



rcg + max{uii } + S — 1 
n- 1 



(2) 



On the other hand, the estimate on the degree of a quasi- 
homogeneous variety gives the following complexity bound 



for FGLM: Cfc 



O n 



, where D — n™=i d% is 



the degree of the ideal {Yiom w (F)) . In the end, for the whole 
process, we can see that the complexity bounds for our direct 
strategy is divided by a P u factor, when compared to the 
strategy of going through hom w . 

A consequence of these results is that in general, the 
higher the weights of the variables, the bigger the gain. In 
particular, it shows why the applying matrix-Fs on hom VK (_Fj 
did not yield satisfactory complexity results. 

On the other hand, it also shows that when possible, re- 
placing powers of the variables with higher weight variables 
allows for faster computations of a Grobner basis. Of course, 
on a generic homogeneous or quasi-homogeneous system, it 
won't be possible, but since in practice, the systems we con- 
sider are not generic, situations may appear where this idea 
proves interesting. We will see an example of that in sec- 
tion 4. 

3.4 Precise analysis of matrix-F 5 

We now fix as goal to follow in a more precise way the 
computations occurring in the matrix-Fs algorithm, and to 
obtain better complexity bounds. For this purpose, we take 
on the computations made in [2, ch. 3], without proving 
them whenever the proof is an exact transcription of the 
homogeneous case. 

Let W = (wi, . . . , w„) be a weight system, and fi, . . . , f m 
a system of quasi-homogeneous polynomials in K[Xi , . . . , X n ] | 
which we assume satisfies the hypotheses HI and H2. We 



denote by (dx, . . . , d m ) the respective W-degrees of the poly- 
nomials fi, . . . , fm, and we will assume them to allow the 
existence of such systems. 
We also denote by: 

• At = K[Xi,. . .,Xi], and A = A n ; 

• Si the integer defined in Lemma 9, and S = S n ; 

• Pi = Ilj = l W h and P — Pn\ 

• h = (fi, . . . , fi), and I = I m ; 

• li = {fi, ■ ■ ■ , ft), and I = I„ 
. Ph = deg(/0 = nj = i (dj/t 

• Di = deg(7i) = ]T_i di 



L m • 

/Wi): 



lj=l 



• dreg the degree of regularity of L (or of h) ; 

• Gi the W-GRevLex Grobner basis of T as given by 
matrix-Fs. 

With these notations, we are going to prove the following 
theorem: 

Theorem 10. Let W — (wi, . . . , w n ) be a weight system, 
and fx, ■ ■ ■ ,f m (m <n) a system of W -homogeneous polyno- 
mials satisfying HI and H2. Then the complexity of quasi- 
homogeneous matrix-Fs algorithm (algorithm 1) is: 



Cf 5 = O - D l - 2 )M d 



The goal is to compute precisely how many lines are re- 
duced in a run of the matrix-Fs algorithm, that is, the num- 
ber of polynomials in the returned Grobner basis. This is 
done by the following proposition, which is a weak variant 
of [3, th. 10]: 

Proposition 11. Let W be a system of weights, and (fx,., 
be a W -homogeneous system satisfying the hypotheses HI 
and H2. Let Gi be a reduced Grobner basis of (fx, ... , fi) 
for the VK-GRevLex monomial ordering, for 1 < i < m. 
Then the number of polynomials of W -degree d in Gi whose 
leading term does not belong to LT(Gi_i) is bounded by b d ,i, 
defined by the generating series 



Bi(z) = ^2 b d> i 



1 



PROOF. The proof of [3, th. 10] still holds in the quasi- 
homogeneous case, using formula (1) for the Hilbert series 
of a quasi-homogeneous regular sequence. □ 

So we can obtain a better bound for the number of ele- 
mentary operations performed in a matrix-Fs run. Indeed, 
Bi(l) represents the number of reduced polynomials in the 
computation of a Grobner basis of (fx, . . . , fi, Xi+x, ■ ■ ■ , X n ), 
that is as many as in the computation of a Grobner basis 
of (fx, . . . , fi): since we only perform reductions under the 
pivot line, [3, prop. 9] shows that the lines coming from 
Xi+i, . . . ,X n will not add any reduction. Note that the 
above generating series is the same as the Hilbert series of 
(fx, ■ ■ . , fi-i, Xi, . . . , X n ), and so, that its value at z = 1 is 
the degree of that ideal, or Di-i. Therefore, we know that 
the number of reduced polynomials with label (m, fi) will 
be Di-i — Di-2 (with convention that Do = 0). 

Now, let g be any polynomial of W-degree d being reduced 
in a run of the matrix-Fs algorithm on (fx, ■ ■ ■ , fi)- From [3, 
prop. 9], we know that the leading term of g, after reduction, 
is in Ai. So overall, in W-degree d, we reduce by at most as 
many lines as there are monomials in Ai, that is Md,w(i)- 
Furthermore, each reduction costs at most 0(Md,w(n)) ele- 
mentary algebraic operations, since this is the length of the 
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Number of operations 




24 30 36 4248 



Matrix-Fs run on hom w (F) 
■ Number of operations 
1 — B — 1 Bound (4) 

- b- Bound (2) with W = 1 

Weighted variant of matrix-Fs (Algo. 1) 

• Number of operations 
— • — Bardet-Faugere-Salvy bound: (5) 
-e- Bound (3) 

- «- Bound (2) 



Figure 1: Bounds and values, on a log- log scale, for the number of arithmetic operations performed in matrix-Fs for a generic 
system with W = (1, 2, 3) and D = (d, d, d) 



matrix lines. And we perform these reductions up to degree 
dreg. Note that, if i — 1, there clearly isn't any reduction in 
the computation, and we obtain the following formulas: 



C Fs = O 



= o 



»=a 

m 

E 



(n) (3) 



A- 



-Pi Pri 



M 4% +Si - 



(i) -M M (n) 

-1,1 v y (4ei+Sn-n+l,l V ' 



In comparison, the above reasoning for matrix-Fs applied 
to F would give 



C Fs = O 



so that here again, working with quasi-homogeneous polyno- 
mials yields a gain or roughly P 3 . Note that the exponent 3 
(instead of the previous w) is not really meaningful, because 
we assumed here that we were using the naive pivot algo- 
rithm to perform the Gauss reduction. However, if we as- 
sume u) = 3 in the previous computations as well, we observe 
that our new bound is generally much better than the pre- 
vious one: figure 1 shows a plot of data obtained both with 
algorithm 1 and with matrix-Fs through hom ff , together 
with the different bounds we can compute. 

Asymptotically, though, the gain does not look important, 
since the complexity is still 0(nD 3 ) where D is the degree of 
the ideal and n > m the number of variables, or in 0(nd 3n ) 
where d is the greatest di. 

Remark 3. One may also push the computations a bit fur- 
ther, and obtain an even more accurate bound, expressed 
in terms of the bd,i (these calculations are done in [2] for 
the homogeneous case, and can easily be transposed to the 
quasi- homogeneous case): 



P l + lPn 



■M d 



• M, 



(n) . (5) 



d+d 



As an example, we computed that bound as well for a par- 
ticular case, and included it in figure 1. As can be seen, this 
bound is indeed better than the intermediate evaluation (3), 
but the difference is low enough to justify using the latter 
evaluation. Furthermore, the bound (3) expressed in terms 
of the Di's is more useful in practice, since it has a closed 



formula using only the parameters of the system (n, m, di 
and Wi). This allows us to use it in complexity evaluations, 
in both theory and practice. 

Remark 4- As can be seen on figure 1, the number of 
operations needed by Matrix-Fs on the homogenized sys- 
tem is not significantly higher than the number of opera- 
tions needed by the quasi-homogeneous variant of Matrix-Fs. 
That is mostly true because the unmodified algorithm can 
make use of some of the structure of the quasi-homogeneous 
systems (for example, columns of zeroes in the matrices). 

4. THE AFFINE CASE 

We will now consider the case of non-necessarily-quasi- 
homogeneous polynomials. One of the methods to find a 
GRevLex Grobner basis of such a system is to apply F5, 
considering at W-degree d the set of monomials having W- 
degree lower than or equal to d. This is equivalent to homo- 
geneizing the system, that is to adding a variable Xi > • • • > X n > ff( 
and applying the classical F5 algorithm to this homogeneous 
system. The reverse transformation is done by evaluating 
each polynomial at H = 1. 

However, this process makes computing the complexity 
of the F5 algorithm harder. The main reason is that deho- 
mogenizing does not necessarily preserve W-degree, and as 
a consequence, it is no longer true that running the matrix- 
F5 algorithm up to W-degree d provides us with a basis, 
truncated at Vl^-degree d. What remains true though is that 
past some W^-degree, we may obtain a Grobner basis for the 
entire ideal. 

To ease the study, we will only consider regular in the 
affine sense systems (as found in [ ']), that is, systems for 
which no such degree fall may happen. 

Definition 2. Let W be a weight system, and (/1, . . . , f n ) 
be a system of non-necessarily W-homogeneous polynomials. 
We write, for any 1 < i < n, hi the quasi-homogeneous com- 
ponent of highest W-degree in /j. We say that the sequence 
(fi) is regular in the affine sense when the sequence (hi) is 
regular (in the quasi-homogeneous sense). We call degree 
of regularity of the ideal (fi) the degree of regularity of the 
ideal (hi). 

Since a degree fall in a run of matrix-Fs is precisely a re- 
duction to zero in the highest W-degree quasi-homogeneous 
components of the system, we know that the F5 criterion 
rules out all degree falls in a run of matrix-Fs on such a 
regular system. 



System 



Generic n = 7, W = 


(1,1, 1, 1,2,2,2), D = (4, . . . ,4) 


2048 


2.7s 


1.2 


0.43s 


2.6 


Generic n = 8, W = 


(1,1,1,1,2,2,2,2), D = 


(4,..., 4) 


4096 


12.3s 


1.8 


2.41s 


3.0 


Generic n = 9, IV = 


(1, 1, 1, 1, 1, 2, 2, 2, 2), D 


= (4,..., 4) 


16384 


314.9s 


2.5 


119.6s 


2.7 


Generic n = 7, W = 


(2,. . . ,2, 1, 1), D = (4, . 


■ A) 


512 


0.09s 


3.2 


0.06s 


1.7 


Generic n = 8, W = 


(2, ... ,2, 1, 1), D = (4, . 


■A) 


1024 


0.39s 


4.2 


0.17s 


1.9 


Generic n = 9, W = 


(2, ...,2,1,1), D = (4, . 


■ A) 


2048 


1.63s 


1.9 


0.59s 


2.0 


Generic n = 10, IV = 


= (2, ...,2,1,1), D = (4, 


■ ■■A) 


4096 


7.54s 


5.4 


2.36s 


2.6 


Generic n = 11, W = 


= (2, ...,2,1,1), D = (4, 


■ ■■A) 


8192 


33.3s 


6.4 


17.5s 


2.4 


Generic n = 12, IV = 


= (2, ...,2,1,1), D = (4, 


■ ■■A) 


16384 


167.9s 


6.8 


115.8s 




Generic n = 13, W - 


= (2, ...,2,1,1), D = (4, 


■ ■■A) 


32768 


796.7s 


8.4 


782.74s 




Generic n = 14, W = 


= (2, ...,2,1,1), D = (4, 


■ ■■A) 


65536 


5040.1s 


CO 


5602.27s 




DLP Edwards n = 4 






512 


0.03s 


3.7 


0.07s 




DLP Edwards n = 5 






65536 


935.39s 


6.9 


2164.38s 


3.2 



Table 1: Benchmarks with FGb on some affine systems 



We can study the complexity of F5 by looking at a run of 
matrix-Fs on the homogenized system. As an example, we 
prove the following theorem: 

Theorem 12. Let W = (wi, . . . ,w n ) be a weight system, 
and fi , . . . , f m a generic system of polynomials of the form 
fi = <?i + Ai, with gi W -homogeneous of W-degre di and 
Xi 6 K. Let D be the degree of the system, d rcg its degree 
of regularity, and 5 the gcd of the di 's. We can compute a 
IV-GRevLex Grobner basis of this system in time 

O ( ^M d , w {nf 



For each system, we ran the F5 and FGLM algorithms on 
both the homogenized system obtained with hom w (£f 5 (qh) 
and tFGLM(qh)) and the unchanged system (for which we 
provide the ratio of computing times). 



or in other words, we can divide the known complexity of the 
F$ process on such a system by 8 U . 

Proof. The idea is that when we homogenize the system, 
we can choose any suitable weight for H, not necessarily 1. 
More precisely, we can set the weight of H to be 8, so that 
the homogenized polynomials become //* ~ gi + XiH di ^ s . 

This way, assuming the computations made at section 2.2 
still hold, we have the same improvements on the bound on 
dreg and on the size of matrices as before, and we thus have 
the wanted result. 

Note that even if the initial system is generic, the homog- 
enized system is not. However, one can check that if the 
initial system was regular in the affine sense, the homoge- 
nized system is still regular. Indeed, it's enough to check 
that no reduction to zero occur in a matrix-Fs run, but it is 
clear, since such a reduction would in particular be a degree 
fall. Also, the property of being in Noether position for the 
m first variables is clearly kept upon homogenizing. 

As such, generically, our homogenized system is regular 
and in Noether position, so the previous computations in- 
deed still hold. □ 

5. EXPERIMENTAL RESULTS 

We have run some benchmarks 1 , using the FGb library [1 1] [ 
We present these results in Table 1. The examples are cho- 
sen with increasing n (number of variables and polynomials), 
two different classes of systems of weights W and systems 
of W^-degrees D. With these conditions, we built a generic 
system of polynomials fi in F6552i[X], such that all mono- 
mials appearing in fi have VF-degree at most di. The last 
examples are systems arising in the study of the Discrete 
Logarithm Problem, when trying to compute the decompo- 
sitions of points on an elliptic curve (see [12]). 

1 A11 the systems we used arc available online on 
http : //www-polsys . Iip6 . f r/~ jcf /Sof tware/benchsqhomog.html J 
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